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Multidrug resistance is a serious threat to public health. Proton motive force-driven antiporters from the major 
facilitator superfamily (MFS) constitute a major group of multidrug-resistance transporters. Currently, no reports 
on crystal structures of MFS antiporters in complex with their substrates exist. The E. coli MdfA transporter is a 
well-studied model system for biochemical analyses of multidrug-resistance MFS antiporters. Here, we report three 
crystal structures of MdfA-ligand complexes at resolutions up to 2.0 A, all in the inward-facing conformation. The 
substrate-binding site sits proximal to the conserved acidic residue, D34. Our mutagenesis studies support the struc- 
tural observations of the substrate-binding mode and the notion that D34 responds to substrate binding by adjusting 
its protonation status. Taken together, our data unveil the substrate-binding mode of MFS antiporters and suggest a 
mechanism of transport via this group of transporters. 
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Introduction 


Drug resistance has become one of the most serious 
threats to public health, with devastating consequences 
for the treatment of both infectious diseases and cancers 
[1-3]. Both eukaryotic and prokaryotic cells have de- 
veloped various evasion mechanisms to overcome the 
toxic effects of chemically distinct drugs. One simple 
and effective way is to employ multidrug-resistance 
transporters with broad substrate specificity to pump 
the drugs into the extracellular space [4]. Many of these 
multidrug-resistance transporters belong to the major 
facilitator superfamily (MFS) [5, 6]. As one of the largest 
secondary active transporter superfamily, MFS trans- 
porters are of broad biological importance, with some of 
them possessing poly-specific substrate-binding ability. 
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They utilize a secondary electrochemical potential to 
transport substrates, often against the concentration gra- 
dient of the latter. MFS transporters are further classified 
as symporters, antiporters and uniporters. Bacterial MFS 
antiporters often utilize proton motive force (PMF), 
including a negative-inside Apy and an alkaline-inside 
Alyx, to drive the expulsion of hydrophobic antibiotics 
[7]. In general, to enter the target cells, drugs usually 
carry positive charges, thus facilitating penetration of 
the cell membrane bearing a negative-inside membrane 
potential. More than often, these positive charge-carrying 
drugs can be expelled by cellular transporters that belong 
to the cation-proton antiporters. Therefore, understand- 
ing the mechanisms of MFS antiporters is of tremendous 
theoretical and practical importance. 

MEFS transporters contain a 12 transmembrane 
(TM)-helix core consisting of two pseudo-symmetri- 
cal, six-helix domains [8, 9]. A membrane-embedded 
central cavity is present between these two domains, 
thus forming the substrate-transport path. In previously 
reported crystal structures, two major conformations of 
MFS transporters have been observed, i.e., inward-facing 
(Cn) and outward-facing (Cow) conformations [8-10]. 


Conformational changes from Cou to Cn and from Cp to 
Cou are thought to utilize a rocker-switch mechanism for 
substrate transport [8]. One of the fundamental questions 
about MFS transporters is the mechanism by which the 
electrochemical potential alteration is coupled with the 
conformational changes. 

Being one of the best-characterized bacterial MFS 
antiporters biochemically, MdfA from Æ. coli (ecMdfA) 
is known to confer resistance to a variety of structurally 
distinct cationic and zwitterionic lipophilic compounds, 
as well as to a number of electroneutral antibiotics of 
clinical importance [11]. Overexpression of MdfA from 
plasmids has been observed in multidrug-resistance E. 
coli strains isolated from clinical patients [12]. MdfA or- 
thologs are found in several pathogenic bacterial species 
(Supplementary information, Figure S1A). Because of 
its broad substrate specificity, MdfA is predicted to have 
a large hydrophobic substrate-binding pocket, which is 
able to bind multiple substrates either individually or 
even simultaneously. In antiporters, substrate binding and 
protonation are often found to compete with each other 
[13], in contrast to symporters which exhibit a positive 
coupling between substrate binding and protonation [14, 
15]. However, fundamental questions remain regarding 
the structures of substrate-binding site(s) and protonation 
site(s), and the functional coupling between these two 
events. 

Four motifs, A-D, have been identified in MFS anti- 
porters [7] (Supplementary information, Figure S1). Two 
of them, motifs A and B, are common in antiporters as 
well as symporters. Motif-C is specific for antiporters 
including MdfA, and is therefore called the antiporter 
motif [16]. Motif-D is only found in some subgroups of 
MFS transporters, but is conserved in MdfA orthologs. 
The function of motif-A in stabilizing the Cow state has 
been elucidated in a recent study reporting the crystal 
structure of a putative antiporter, Æ. coli YajR [17]. Thus, 
motif-A is unlikely to contribute to substrate recognition. 
Motif-B was proposed to be involved in energy cou- 
pling [18] and substrate binding-induced conformational 
change [19], yet the precise mechanisms for its actions 
remain elusive. In MdfA, motif-D contains two acidic 
residues essential for transport activity [13, 20]. Current- 
ly, structural studies on PMF-driven MFS antiporters are 
limited. To date, no crystal structures have been reported 
for MFS antiporters in complex with their substrate(s). 
The antiporters with known crystal structures, EmrD and 
YajR [17, 21], lack detailed biochemical characterization 
[22]. Furthermore, the functional roles of motifs B, C 
and D in the energy coupling processes of multidrug-re- 
sistance MFS antiporters remain to be elucidated. 

Here, we report the crystal structures of ecMdfA in 
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complex with a known substrate chloramphenicol (Cm) 
as well as substrate analogs that are structurally distinct, 
i.e., deoxycholate (Dxc) and n-dodecyl-N,N-dimethyl- 
amine-N-oxide (LDAO), all in the C,, state. As an antibi- 
otic, Cm binds to the 50S subunit of bacterial ribosome 
and inhibits peptide synthesis through interactions with 
nucleotides of the rRNA [23]. All four motifs of MFS 
antiporters are illustrated in the three-dimensional (3D) 
structure of ecMdfA. Based on the structure, we verified 
the functional roles of motif-D, including the two con- 
served proton-titratable residues inside the central cavity, 
E26 and D34. Our analysis provides structural evidence 
that D34 is the acidic residue whose deprotonation in 
the C,, state is more affected by both motif-B and sub- 
strate binding than E26. Potential functions of motif-B 
in regulating the deprotonation process and the competi- 
tion between substrate binding and protonation are also 
discussed. Furthermore, isothermal titration calorimetry 
(ITC) methods were used to confirm the Cm binding 
mode. 


Results 


Overall structure 

The £E. coli MdfA consists of 410 amino acid residues. 
Three crystal structures of MdfA were determined in 
complex with different ligands, i.e., Cm (2.4 A), Dxc 
(2.0 A), and LDAO (2.2 A), all of which exhibited MdfA 
in the C,, state (Supplementary information, Table S1). 
Among these ligands observed, Dxc was originally 
used as an additive for successful crystallization, and 
LDAO was used as a component of a detergent mixture 
employed for protein purification and crystallization. 
The MdfA-Cm crystal was obtained by soaking the 
MdfA-Dxc crystal in a solution (pH 5.8) supplemented 
with 5 mM Cm. The MdfA-LDAO crystal was obtained 
when the pH of the crystallization conditions was raised 
from 5.8 to 8.5. In all three crystal structures, the final 
refined model contained the intact peptide chain of res- 
idues 9-400. Except for the ligands, the three structures 
are nearly identical, with root mean square deviations 
ranging from 0.3 to 0.6 A (for all available Ca atoms). 
We therefore focus on the structure of the MdfA-Cm 
complex as a representative of all structures, unless spec- 
ified otherwise. Overall, the crystal structure of MdfA 
consists of 12 TM helices (TMs 1-12), with TMs 1-6 in 
the N-terminal part and TMs 7-12 in the C-terminal part 
(hereinafter referred to as N-domain and C-domain, re- 
spectively), and a 12-residue, inter-domain, amphipathic 
helix (04.43 Figure 1). The MdfA structure consists of 
four pseudo-symmetrical, three-helix repeats, consistent 
with the typical structure of an MFS transporter. The 
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Figure 1 Overall structure of MdfA-Cm. The backbone of MdfA is shown in tube representation, with N-domain, C-domain 
and as in green, cyan and orange, respectively. The substrate Cm is shown as magenta stick model. Positions of E26 and 
D34 are marked as red spheres, and R112 as a blue sphere. The inward-facing cavity is shown as dot-surface representa- 


tion. TMs are labelled in the right panel. 


first helix in each repeat (TMs 1, 4, 7 and 10) contrib- 
utes to the formation of the central cavity, and is thus 
referred here to as the cavity helix. The second helix in 
each repeat (TMs 2, 5, 8 and 11) is long and curved and 
is referred to here as the rocker helix. The third group of 
helices (TMs 3, 6, 9 and 12) are located at both ends of 
the longest dimension of the TM core, referred to here as 
support helices (Figure 1; a more detailed description of 
the overall structure can be found in the Supplementary 
information, Data S2). In order to specify the position of 
a given residue relative to the membrane, we divided the 
membrane-embedded part of each TM helix (ca. 20 resi- 
dues or five helix turns) into five regions, and numbered 
them as 1-5, starting from the periplasmic side [24] (Sup- 
plementary information, Figure S2). For example, TM7.3 
refers to the central region of TM 7. Such a description is 
helpful for discussion of the interaction between a mem- 
brane protein and its membrane environment. 


Titratable residues inside the cavity 

There are only two proton-titratable residues inside the 
central cavity of ecMdfA, E26™'* and D34™'? (Figure 
2A), both located in motif-D. They are the candidates of 
protonation sites and are known to play important roles 
in the transport activity of MdfA [13]. In the Cn state, 
E26 is located on the interior wall of the central cavi- 


ty and is surrounded by Yao", Liga, Milado ns, 
v149™>? and A150'™°°, and its solvent-accessible sur- 
face (SAS) area is 5 A’ (using a 1.4-A radius probe). 
D34 is located at the apex of the cavity, is surrounded 
by N33", M58" A153", P54”? and 12397", 
and is about 10 A (the shortest distance) away from E26. 
Its SAS area is 12 A’ in the C, state. In comparison, 
a completely solvent-exposed aspartate residue, e.g., 
D343, has an SAS of over 120 A’. Thus, these two acidic 
residues are buried mostly by hydrophobic residues from 
the N-domain. In addition, the putative Co,, structure 
was modelled based on the C, structure of MdfA and the 
Cou structure of YajR (PDB ID: 3WDO; Supplementary 
information, Figure S3). In the putative Cow state model, 
the SAS area of D34 increases to about 28 A’ because of 
movement of 1239 of the C-domain, while that of E26 
remains the same as in the C,, state. 


Substrate binding and motif-D 

In the structure of MdfA-Cm, Cm was found inside 
the central cavity which is ~3 000 A’ in size (Figure 2B 
and Supplementary information, Figures S4A and S5). 
Twelve amino acid residues are involved in the binding 
to Cm (Supplementary information, Table S2). These 
surrounding residues come from both cavity helices (TMs 
1, 4 and 7) and rocker helices (TMs 2 and 5), with most 
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Figure 2 Structural and functional studies of MdfA in Cm binding. (A) Stereo view of the titratable residues. E26, D34 and 
their surrounding residues are shown as stick models. Distances between selected atom pairs are marked as dash-lines and 
labelled. Backbone of the protein is shown in wheat tubes. The substrate Cm is removed for clarity. (B) Stereo view of the Cm 
binding. The substrate Cm and its surrounding residues are shown as stick models. H-bonds shorter than 3.0 A are shown as 
blue dash-lines. Electron density of the omit map is contoured at 3 o (blue) and 2 o (wheat) levels. The views of A and B are 


different by ~180°. (C) ITC analysis on Cm binding. 


of them displaying some degree of conservation among 
MdfA orthologs (Supplementary information, Figure 
S1A). Among them, Y30™!'?, N33™'? and D34™!? 
come from motif-D, and the rest are all with hydrophobic 
side chains. Three hydrogen (H) bonds are observed be- 
tween Cm and the side chains of N33 and D34. In partic- 
ular, the carboxyl side chain of D34 forms two H-bonds 
with O4 and O5 hydroxyl groups of Cm (Figure 2). 
Since a hydroxyl group can be either an H-bond donor or 
acceptor, both protonated and deprotonated forms of D34 
may maintain the H-bonds. In contrast, E26™'* is not 
directly involved in Cm binding. In addition, the nitryl 
group of Cm points to the solvent space of the cavity. It 
may explain the observation that thiamphenicol, which 
has a substitution at this nitryl group, is also a substrate 
of MdfA [5] (Supplementary information, Figure S5). 
Besides Cm, we also solved the crystal structures 
of two other MdfA-ligand complexes, MdfA-Dxc and 
MdfA-LDAO (Supplementary information, Figures 
S4B, S4C and S5). Dxc binding was shown to stabilize 
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ecMdfA in an in vitro thermofluor assay (Supplementary 
information, Figure S6). As mentioned above, LDAO 
was found in the central cavity when the pH of the 
crystallization conditions was adjusted from 5.8 to 8.5. 
Deprotonation of D34 at higher pH may be responsible 
for the binding of the positively charged head group of 
LDAO. A substrate binding pattern similar to that of Cm 
was observed in the crystal structures of the Dxc and 
LDAO complexes (Supplementary information, Fig- 
ure S4D and Table S2). In particular, residues Y30"™'”, 
N33™!?, D34™!? and L236'™”” are involved in ligand 
binding in all three crystal structures. The substrates sit 
between D34™!? and L236™”°, and have more contacts 
with the N-domain than with the C-domain (Figure 2B). 
Reduction of the solvent-accessible volume of the cavity 
varies upon the binding of different ligands: Cm, ~450 
A? (equivalent to ca. 15 water molecules); Dxc, ~700 
A’ (23 H,O); LDAO, ~490 Å? (16 H,O). The exclusion 
of solvent molecules from the central cavity would de- 
crease the dielectric constant, enhance the electrostatic 
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field inside the cavity, and affect the protonation status of 
titratable residue(s). In addition, in the LDAO complex, 
the positively charged group of LDAO is located close to 
D34, while in the Dxc complex the potentially negatively 
charged group of Dxc is distal from D34. Therefore, our 
structural analyses of MdfA-ligand complexes suggest 
that MdfA may accommodate chemically distinct sub- 
strates. 

We next used the ITC method to test the mutagenesis 
effects of Y30, N33 and D34 on Cm binding (Figure 2C). 
WT MdfA showed a K, of 75 (+ 8) uM for Cm, while 
Y30A, N33A and D34A mutants lost most of the binding 
ability towards Cm. A previous report showed that D34X 
variants maintained TPP’ binding ability while losing the 
transport activity [13]. The binding mode of TPP’ is like- 
ly to be different from that of Cm, in which TPP’ may 
not directly bind with D34. Moreover, TPP’ stabilizes 
MdfA in its Co, state in the in vitro binding assays, dis- 
tinct from the C, state observed in our crystal structures. 
Taken together, these binding analyses strongly support 
our structurally observed Cm-binding mode. 


Motif-B in the 3D structure 

Motif-B (i.e., “R'’’xxQG”, where capital letters stand 
for highly conserved residues) and its surrounding resi- 
dues in the 3D structure are conserved among MdfA or- 
thologs (Figure 3 and Supplementary information, Figure 
S1A). In particular, the essential basic residue R112 par- 
ticipates in an H-bond network which includes C96™*”, 


Periplasm 7 


N-domain <“ 


E26 2 ™ 


Cytosol 


Qi1s™*?, R112™*?, G32™!? (main chain), N33™’* and 
a well-ordered water molecule. This H-bond network 
is further insulated by other hydrophobic residues, e.g., 
F174™®** (g-orbit) and Y61'* (x-orbit). Moreover, the 
distance between R112 and D34 is 9 A, and that between 
R112 and E26 is 16 A. Thus, motif-B is likely to affect 
D34 more than E26. Together, the “R'’xxQG” motif in 
TM4 and conserved surrounding residues will be referred 
to further on as “3D motif-B”. As will be discussed lat- 
er, we postulate that this motif plays an essential role in 
coupling the protonation status with substrate binding in 
MdfA. 


Drug resistance assays 

To examine the functional relevance of a number of 
MdfA mutations, we performed Cm-resistance assays. 
Bacterial cells expressing WT MdfA showed resistance 
to Cm (2.5 pg/ml) in solid medium, while control cells 
containing the empty pET28a vector did not (Figure 4A). 
Expression levels of WT MdfA and all mutant variants 
were verified by immunoblotting (Supplementary in- 
formation, Figure S7A). Mutations at the two titratable 
residues of motif-D, E26A/Q and D34A/N, abolished the 
resistance, as previously reported [25], and Y30A/F mu- 
tants with mutations also in motif-D region partially lost 
the Cm resistance (Figure 4B and 4C). Together, these in 
vivo results support the structural observation that D34 
and Y30 are directly involved in Cm binding. 

Furthermore, to verify the Cm-binding mode observed 
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Figure 3 Stereo view of motif-B in the 3D structure. Residues of motif-B and surrounding conserved residues are shown as 
stick models. Selected Gly residues (Ca atoms, wheat coloured) and water molecules (red) are shown in spheres. Backbones 
of the N- and C-domains are shown in tubes, coloured in wheat and red, respectively. 
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Figure 4 Drug resistance of MdfA variants. For each variant of MdfA, a serial dilution of cell culture was spotted on solid 
medium containing 0.2 mM IPGT and 2.5 ug/ml Cm. Photos were taken 36 h after incubation at 37 °C. Experiments were 
repeated at least three times, and the results were reproducible. Representative results are shown (also see Supplementary 
information, Figure S7). Positions of the mutation sites are labelled in Figures 2A, 2B and 3. 


in the crystal structure, we systematically mutated the 12 
residues that are involved (Figure 2B and Supplemen- 
tary information, Table S2). In particular, non-alanine 
residues were mutated to alanine, and A150 was mutated 
to leucine. Most of these point mutations resulted in a 
partial or complete loss of Cm resistance (Figure 4C), 
supporting the structurally observed Cm-binding mode. 
In addition, we analyzed the effects of mutations in 
the 3D motif-B on Cm resistance. We first examined 
mutations at R112 and Q115, which form an H-bond in 
the MdfA-Cm structure. The point mutations R112Q and 
Q115E/N/R indeed resulted in loss of Cm resistance, 
whereas the R112H strain (which maintains the positive 
charge at position 112) was able to retain partial resis- 
tance (Figure 4D). Furthermore, the double mutation 
R112Q/Q115R, which switches the two interacting res- 
idues R112™** and Q115'™**, did not bestow Cm-resis- 
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tance. In addition, an extra mutation in close proximity 
to R112, either G32R or S57R, was introduced together 
with R112C to investigate whether the previously re- 
ported inactive mutation R112C could be rescued [26]. 
G32R/R112C indeed restored Cm resistance, while 
S57R/R112C failed, probably due to improper position- 
ing of the critical positive charge. To analyze residues 
surrounding motif-B, mutations C96A/T/S were intro- 
duced. The result showed that C96S was functional, 
while C96A/T abolished Cm resistance (Supplementary 
information, Figure S7C). The conserved residues Y61, 
F108 and F174, proposed above as part of the 3D mo- 
tif-B (Figure 3), were mutated. Among these mutants, 
Y61L, F108L/V and F174L lost Cm resistance, and only 
Y61F partially maintained the activity (Supplementary 
information, Figure S7C). N33, which is involved in Cm 
binding, is located between the key residue of motif-B, 
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R112 and the proton-titratable residue D34. The N33L 
mutant exhibited a loss of Cm resistance, whereas N33D 
maintained the activity (Figure 4B). Together, these 
results suggest that the function of MdfA requires main- 
taining the integrity of the 3D motif-B, including the 
presence of the essential positive charge (R112), the bur- 
ied H-bond network, as well as the surrounding aromatic 
residues. 

More mutational analyses were also performed to 
investigate the functional roles of motifs A and C and to 
confirm our Cou model of MdfA (Supplementary infor- 
mation, Figure S3B). These results are shown in Supple- 
mentary information, Figure S7. 


Conformational changes of MdfA upon substrate binding 

In MFS antiporters, a proline cluster in TM5 (i.e., mo- 
tif-C of MdfA) is considered to be a signature motif [27]. 
In the C, crystal structure of MdfA, motif-C is found to 
be located in the interface between the N- and C-domains 
(hereinafter referred to as the inter-domain interface) and 
close to the substrate binding site. To probe the confor- 
mational changes upon substrate binding and to elucidate 
the role of motif-C in the alternating access mechanism 
of the H'-drug antiporter, we used the smFRET imaging 
method [28] to investigate conformational changes of 
MdfA as well as its motif-C mutants. 

First, we constructed a Cys-less derivative of MdfA 
(denoted as WT*) by mutating all four Cys residues in 
MdfA to either Ser or Ala, and confirmed the WT-like 
function of WT*-MdfA in the Cm-resistance assay (Fig- 
ure 5A). Then, multiple cysteine pairs were individually 
introduced into the N- and C-domains as potential la- 
belling sites of fluorescent dyes, and the Cm-resistance 
activities of these MdfA variants were verified. Among 
them, WT*-MdfA carrying the additional V167C/N374C 
Cys-pair mutations (~33-A apart) on the periplasmic 
side was chosen (denoted as WT**), due to its better 
signal-to-noise ratio, to probe the conformational status 
of MdfA in response to Cm binding. In particular, the 
detergent-solubilized, dye-labelled WT**-MdfA sample 
showed a higher population of the C, state than the Cow 
state at pH 8.0 (Figure 5B, where a lower FRET signal 
indicates a longer distance between the paired probes). 
Furthermore, effects of drug binding on the conforma- 
tional change were investigated. The result showed that 
Cm binding induced closure of the periplasmic side of 
MdfaA. This result is in agreement with a previous report 
showing that, in the absence of electrochemical poten- 
tial, binding of Cm and EtBr stabilizes the C,,, state while 
TPP’ binding stabilizes the C,,, state [29]. 

In addition, we investigated the effects of point mu- 
tations in motif-C on conformational populations. The 


results showed that both P154A and P158A variants were 
prone to an outside-open (Cow) conformation in the ab- 
sence of Cm. Similar to WT**, addition of Cm shifted 
the MdfA conformation to an outside-closed (Cp) state 
(albeit to a less extent than WT**). Based on the C,, 
crystal structure, we reasoned that these point mutations 
induced conformational changes in the inter-domain in- 
terface, thus destabilizing the C, state and favouring the 
Coy State. 


Discussion 


Multidrug resistance 

MdfA is a paradigm of multidrug-resistance MFS 
antiporters [11]. Our structural analysis at resolution 
range of 2.0-2.4 A provides direct information on its 
substrate-binding mode in the C,, state. Of about 30 res- 
idues that form the substrate-binding cavity in the C,, 
state, over two-thirds are hydrophobic. This observation 
is consistent with the fact that most substrates of MdfA 
are hydrophobic in nature [11]. Further, using their polar 
(or positively charged) groups, all three ligands shown 
in the crystal structures bind to the apex of the cavity, 
which is located close to the titratable residue D34 (Fig- 
ure 2). Their hydrophobic parts pack against the hydro- 
phobic cavity wall from the N-domain as well as the 
inter-domain interface between TMs 2 and 11 (Figure 1). 
Moreover, among these three ligands, Cm is a bona fide 
substrate of MdfA, and the other two share recognizable 
structural features with a number of known substrates of 
MdfA [5] (Supplementary information, Figure S5). In- 
terestingly, it was found recently that Dxc is a substrate 
of another E. coli MFS transporter, MdtM [30], which 
shares 41% sequence identity and 63% similarity with 
MdfA (Supplementary information, Figure S1B). In 
our Dxc-MdfA crystal structure, we found that residues 
involved in Dxc binding are conserved between MdfA 
and MdtM. Thus, Dxc is a potential substrate of MdfA. 
The presence of these ligands in the complex crystal 
structures likely reflects the in vivo substrate preferences 
of MdfA. Moreover, ligand-binding modes of poly-spe- 
cific multidrug-sensing proteins of the soluble TetR re- 
pressor family have been well studied, e.g., in the Cm- 
TtgR complex crystal structure [31]. Although detailed 
ligand-binding modes of TetR proteins and MdfA are 
different, a common feature of the ligand-binding cav- 
ities is that hydrophobic resides line up the side walls 
of the cavity. In the cases of the TetR family members, 
the common function of ligand binding is to induce a 
conformational change of the protein homodimer, thus 
triggering the disruption of these repressor dimers. Al- 
though we showed the detailed binding modes of MdfA 
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with different ligands, we would like to argue that, for a 
multidrug-resistance antiporter, the change of the proton- 
ation status of the titratable residue (e.g., D34 in MdfA) 
is more critical for the conformational change of the 
transporter (thus the transport activity) than the detailed 
substrate-binding modes. More specifically, the interac- 
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Figure 5 smFRET data. (A) Cm resistance of WT* and WT**. 
C43 cells expressing the two MdfA variants, Cys-free back- 
ground (WT*) and V167C/N374C/WT* (WT**), showed similar 
level of drug resistance to 2.5 ug/ml of Cm as the WT. (B) FRET 
histograms from single-molecule traces of WT**, P154A/WT** 
and P158A/WT**. Note that protein samples used in this ex- 
periment were incubated in DDM (0.05%), pH 8.0. Data were 
collected in the absence (left) and presence (right) of Cm (100 
uM). Results of two-peak, nonlinear Gaussian fit (with the soft- 
ware Origin) are shown as solid lines in cyan and red, repre- 
senting the outside-open (Cou) and outside-close (Cp) states, 
respectively. The black dotted line is for the summation of the 
two Gaussian distributions. 
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tion between the protonation site and membrane potential 
plays an important role in determining whether the C,, or 
Cou State is thermodynamically favoured [17]. 

Competition between substrate binding and protonation 
is likely to be a universal mechanism for all PMF-driven 
antiporters, although it may vary in detailed implemen- 
tation among different antiporter families [32, 33]. In the 
central cavity of MdfA, there are only two proton-titrat- 
able amino acid residues, E26 and D34. Both of them 
have been shown to be essential for the transport activity 
[13]. Our crystal structure of Cm-MdfA illustrates that 
D34 is directly involved in Cm binding, while E26 is 
over 8 A away from the substrate (Figure 2). Important- 
ly, the protonation status of D34 may have little effect on 
Cm binding (Figure 2B), yet substrate binding may have 
strong influence on the protonation status (see below). 
Thus, D34 is likely to be the proton-titratable residue 
in MdfA that gets deprotonated upon substrate binding, 
in consistence with the previous finding [13]. Never- 
theless, existence of acidic residues per se inside the 
central cavity, rather than their exact positions, has been 
shown to be essential for transport activity of MdfA [34]. 
For example, either of the two single-point mutations 
T29™'°E and A150'°°E is able to rescue the function 
of the inactive double-mutant variant E26T/D34M [25]. 
V335'™°E is able to re-establish cationic drug resis- 
tance of the inactive mutant E26T [35]. Mutant variants 
L119™® E, V231™*E, G354™"“E and M358™"'°E 
are able to establish divalent cationic drug resistance in 
MdfA, which otherwise lacks such capability due to the 
exchange of only one proton per transport cycle [36]. In 
the Cp structure of MdfA, all residues described above 
except for V231 face the central cavity. V231 is buried, 
but may interact with N331'™'°° which faces the cavity 
(Supplementary information, Figure S8). All these res- 
idues are located at the cytosolic side of the positively 
charged motif-B, and the significance of this distribution 
will be discussed later. 

Substrate specificity of an MFS transporter is achieved 
not only by the matching of shapes between binding site 
and substrate, but also by properly coupling the pro- 
tonation/deprotonation process with substrate binding. 
PMF-driven MFS antiporters appear to be able to medi- 
ate two types of transport, electrogenic and electroneu- 
tral transport. The transport of electroneutral substrates 
(e.g., Cm) by MdfA is electrogenic, while the transport 
of mono-cationic substrates (e.g., TPP’) by MdfA is elec- 
troneutral [37]. In the electrogenic transport, substrate 
binding-induced deprotonation, as suggested by our 
MdfA-Cm crystal structure, would disrupt charge-mem- 
brane potential balance, thus triggering the C,,-to-Co,, 
conformational change. In the cases of mono-cationic 
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Figure 6 Schematic diagram of a mechanistic model depicting the transport cycle of MdfA. 


substrates (e.g., TPP’), substrate binding-induced depro- 
tonation is equally (if not more) important for the Cm- 
to-C,,, conformational change. Without deprotonation, 
export of the positively charged substrate would be 
against the negative-inside membrane potential, which 
is thermodynamically unfavourable. Nevertheless, be- 
cause of the zero net-charge change during the transport 
cycle, membrane potential could not be the direct driving 
force, and other sources of energy are required to drive 
the conformational change. For example, TPP” binding 
has been shown to stabilize the Cow state in addition to 
triggering the deprotonation of D34 [13]. Thus, the dif- 
ferential binding energy between the two states, AGp(s) 
= RT/n(KgoulS)/Kan(S)), is negative, thermodynamically 
favouring the C,,-to-C,y,, conformational change. Taken 
together, a common scheme of PMF-driven antiporters 
is that substrate binding induces deprotonation, while 
different types of substrates may employ distinct mecha- 
nisms coupling the driving force and the conformational 
changes. 


Conserved motifs 

Four motifs, A-D, have been identified in MFS anti- 
porters [7], and all of them are present in MdfA ortho- 
logs (Supplementary information, Figure S1A). Here, 
we focus on the functional relevance of motifs B, C and 
D. Since motif-A functions to stabilize the Cow confor- 
mation [17], motif-A is not in its functional state in our 
MdfA crystal structures, which are in the C,, conforma- 
tion. Discussion on motif-A can be found in the Supple- 
mentary information, Data S2. 


Motif-B According to the PFAM database [38], motif-B 
(“R'’xxQG”) in TM4 is the only high-frequency site that 
contains a membrane-embedded basic residue among 


MFS members (Supplementary information, Figure S2). 
Since the whole 3D motif-B is embedded in the level-2 
region of the N-domain (Figure 3), it is unlikely to be di- 
rectly involved in inter-domain conformational changes. 
Previous mutagenesis analysis of MdfA revealed that, 
among mutants carrying point mutations at position 112, 
only the R112H mutant retained resistance to most sub- 
strate drugs, while R112K maintained some activity and 
R112M/Q/E lost resistance to all antibiotics tested [39]. 
The results of our Cm-resistance assays for mutations in 
the 3D motif-B region (Figure 4D and Supplementary 
information, Figure S7C) agree with the previous report. 
In the putative Cow state, levels 1 and 2 of both cav- 
ity-helices (i.e., TMs 1, 4, 7 and 10) and rocker-helices 
(i.e., TMs 2, 5, 8 and 11) are accessible to solvent from 
the periplasmic side (Supplementary information, Figure 
S3B), thus the positive electrostatic field of motif-B is at- 
tenuated by the solvent which is characterized by a high 
dielectric constant. In the C, state, in contrast, levels 1 
and 2 of the TM helices become buried by the inter-do- 
main interface, and thus the positive electrostatic field 
becomes stronger relative to the Cou state. This enhanced 
electrostatic field is able to promote deprotonation in the 
central cavity, e.g., at D34™'? of ecMdfA. Furthermore, 
in antiporters, loading of hydrophobic substrates from 
the cytosolic side (or directly from the lipid bilayer) may 
strengthen the electrostatic field by reducing the dielec- 
tric constant inside the cavity and/or by moving a posi- 
tive charge into the cavity, thus promoting deprotonation 
(Figure 6). Our structural observations of the substrate 
complexes, particularly the binding of the substrate in 
the vicinity of D34 (Figure 2), are consistent with such 
a mechanism. Taken together, the positive electrostatic 
field generated by motif-B nearby titratable residues in 
the C,, state seems essential for the transport activity of 
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motif-B-containing antiporters. 


Motif-C Motif-C in the rocker-helix TM5 (“GX,GX,G- 
PX,GPX,G”) is also called the antiporter motif [27]. 
In MdfA orthologs, motif-C has evolved into the con- 
sensus sequence “itALM ANvaL iaP'“LI GP’*LvG” 
(Supplementary information, Figure S1). The C-terminal 
(periplasmic) half of motif-C is more conserved than the 
N-terminal half (Supplementary information, Figure S1) 
and packs against a conserved region in TM1 (near P38- 
G39™"" of motif-D). In addition, in the C, crystal struc- 
ture of ecMdfA, L151, P154, L155, P158, and L159 of 
motif-C are directly involved in the inter-domain hydro- 
phobic interaction between TM5 and TM8 (Supplemen- 
tary information, Figure S8). Thus, TMS, in particular 
motif-C, is surrounded by other TM helices, a constella- 
tion that prevents proton leakage in the C, state. 

In motif-C, P158™*' is strictly conserved among 
MdfA orthologs (Supplementary information, Figure 
S1). In general, a proline residue in a TM helix intro- 
duces a break in the helix H-bond network, allowing the 
TM helix to form a kink in the membrane-embedded 
region. The presence of such a proline residue also gen- 
erates a free carbonyl group from the peptide backbone 
of the previous helix turn, allowing the binding of the 
latter with either the side chains of other polar residues 
or embedded water molecules. For examples, in the ec- 
MdfaA crystal structures, the prolyl ring of the conserved 
P154™5? is close to the critical acidic residue D34™'* 
(Figure 2), and P38™'', P158™*! and P243™”" partici- 
pate in the formation of a cage accommodating a cluster 
of four well-ordered water molecules in the inter-domain 
interface (Figure 3 and Supplementary information, Fig- 
ure S8). In agreement with our structural observations, 
single-point mutants P154A and P158A displayed a com- 
plete loss of Cm resistance (Supplementary information, 
Figure S7D). Our smFRET experiment further indicated 
that these mutations shift the conformation population 
to the Cou state (Figure 5B). Similar mutations in the 
TetL antiporter have been shown to cause ion leakage 
[40]. Since C,, is the state in which an antiporter binds 
substrate, it can be argued that it is the default state, and 
thus it is reasonable to suggest that the stability of the 
C,,, State is more important than that of the Cou state for 
an antiporter. Taken together, we propose that motif-C 
has specifically evolved to stabilize the hydrophobic in- 
ter-domain interaction of MFS antiporters in the C,, state. 


Motif-D Among MdfA orthologs, motif-D in TM1 
displays a clear pattern of “E*°fxxY*ianD™migP**g” 
(numbering according to ecMdfA; Supplementary infor- 


mation, Figure S1A). Of the two acidic residues, D34 is 
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absolutely conserved among MdfA transporters, and E26 
is conserved as an acidic residue in about 80% of known 
sequences [13]. Numerous point mutations in motif-D 
have been tested experimentally, showing that this mo- 
tif is important for MdfA functions [13, 20, 25, 35]. It 
has been reported that only one proton is consumed per 
MdfA transport cycle, no matter whether an electro- 
neutral or a mono-cationic substrate is transported [13]. 
However, both E26 and D34 appear to be (de)protonated 
at certain time point. The precise functional roles of these 
two acidic residues seem to be a key to understanding the 
MdfA transport cycle. 

Based on the C, crystal structure of MdfA-Cm (Figure 
2), together with the putative C,,, state model (Supple- 
mentary information, Figure S3) and functional studies 
(Figure 4 and Supplementary information, Figure S7), 
we propose a mechanism in which a proton transfers be- 
tween E26 and D34. Since only one proton is consumed 
per transport cycle, only one residue, either E26 or D34, 
can possibly get protonated at any given time. This sug- 
gests that a proton might be transferred between the two 
residues during the conformational change of MdfA, 
and that the direction of the transfer depends on changes 
of their pK, values. In the Cou state, E26 is at the apex 
of the outward-facing cavity and is less accessible to 
the bulk solvent than D34 (Supplementary information, 
Figure S3B). Thus, E26 is more likely to get the proton 
in the Cou state, although low-probability protonation at 
D34 may occur stochastically. By interacting with the 
negative-inside membrane potential (AV), the protonat- 
ed E26 would provide part of the energy for the Couto- 
C,, conformational change. Loss-of-function mutations 
at E26 can be rescued by a number of point mutations 
in the inter-domain interface of the putative Cou state, 
presumably by decreasing the stability of the latter [29]. 
It suggests that, in WT MdfA, part of the electrostatic 
energy associated with protonated E26 in the presence of 
membrane potential is used to overcome the energy bar- 
rier in the transition state. In the C, state, D34 is located 
at the apex of the inward-facing cavity and becomes less 
accessible to the bulk solvent from the cytosol than E26 
(Figure 2). Thus, in the Cp state D34 becomes protonat- 
ed in the absence of a substrate [13], while E26 becomes 
deprotonated because it becomes closer to the bulk 
solvent. The protonated D34 interacts with the electric 
field of the membrane potential, taking the role of E26 
in maintaining the C, state, and waits for deprotonation 
triggered by substrate binding. Therefore, during the 
Co,,-to-C;, transition, there appears to be a proton transfer 
from E26 to D34; however, it remains unclear whether 
the transfer occurs in a direct or indirect fashion. 

Such a hypothetical proton transfer is in agreement 
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with the following observations: (i) Conserved Y30 
is located between E26 and D34 (Figure 2). It is pos- 
sible that Y30 serves as a proton-wire between these 
two acidic residues. Such a proton-wire could facilitate 
proton transfer in both directions. It might provide an 
explanation to the previous in vitro observation that E26 
facilitates the deprotonation of D34 [13]. (ii) The point 
mutation V335E rescues the loss-of-function mutation 
E26T [35]. Although V335™"°° is located in the C-do- 
main (Supplementary information, Figure S8), it is posi- 
tioned at the apex of the cavity in the putative Cow state 
and becomes solvent-exposed in the Cp state, similar to 
E26 from the N-domain. (iii) The loss-of-function double 
mutation E26T/D34M can be rescued by an additional 
point mutation A150E, and A150'™°° is located at a TM 
level between E26™'* and D34™!* (Figure 2A). In this 
case, one acidic residue at position 150 may play the 
roles of both E26 and D34, i.e., accepting a proton in 
the Cou state, thus driving the Cow-to-C conformational 
change, as well as undergoing substrate binding-triggered 
deprotonation in the C state. Alternatively, the proton 
transfer may not require a proton wire. It is possible 
that the “active” proton moves from the apex of the out- 
ward-facing cavity (e.g., E26 in ecMdfA) to the apex of 
the inward-facing cavity (e.g., D34) simply by dropping 
a proton into and picking another one from the solvent 
pool inside the cavity. 

As another possibility, E26 may facilitate the depro- 
tonation of D34 using the following mechanism. In the 
Cn state, in the absence of a substrate, the electrostatic 
interaction between E26 and D34 is attenuated by the 
solvent inside the cavity. Upon the binding of a hy- 
drophobic substrate, the dielectric constant inside the 
cavity is reduced, and the electrostatic interactions are 
enhanced. Under such a condition, E26 may cooperate 
with R112 and exert an electrostatic attraction from the 
cytosolic side to facilitate the deprotonation of D34. 
Nevertheless, the electronegativity of E26 in the C,, state 
cannot be too strong; otherwise, E26 itself would become 
protonated by attracting a proton from the bulk solvent. 


A putative transport mechanism for MdfA 

Based on available structural and functional informa- 
tion for MdfA, we propose a mechanism for the transport 
cycle (Figure 6). First, in the Cn state and before substrate 
loading, D34 remains protonated (1.e., PKn > pH") 
[13]. This Cn state may be in thermodynamic equilibrium 
with the protonated, substrate-free, Cow state, with the 
Cn state being strongly favoured (such a possible equi- 
librium could explain electroneutral transport via MdfA.) 
At the same time, the C, state is ready for substrate 
loading-triggered deprotonation. By reducing the dielec- 


tric constant inside the cavity, binding of a hydrophobic 
(electroneutral) substrate increases the positive electro- 
static field generated by R112 of motif-B (Figure 6). The 
enhanced positive electrostatic field decreases pK, p34) 
m triggering deprotonation of D34. Upon deprotonation, 
the electrostatic force of the proton-A'Y interaction disap- 
pears, the mechanical balance in the C,, state is broken, 
and some “elastic” energy stored in the C, state is re- 
leased, which then drives the C,,-to-Co,, conformational 
change. Second, in the Cou state, releasing of the sub- 
strate should promote protonation of E26 (Supplementary 
information, Figure $3B). The interaction between the 
protonation site(s) and the negative-inside AY facilitates 
the conformational change of Co,,-to-C,,, and part of this 
input energy is stored as the “elastic” conformational en- 
ergy to be used later [17]. Since protonation provides the 
driving force for the Co,,-to-C,, conformational change, 
substrate releasing is likely to be prior to the protonation 
in order to prevent undesirable backward transfer of the 
substrate. In short, a cycle of electrogenic transport is 
driven by the interaction between protonation and mem- 
brane potential. 

Bibi and co-workers have previously proposed a mod- 
el for MdfA-catalyzed substrate/proton transport based 
on biochemical analyses of proton and TPP” binding 
by MdfA [13]. Although our structural observations are 
mostly consistent with their findings, our newly proposed 
model differs from the previous model [13] in several 
aspects: (i) The previous model considered D34 as the 
protonation site in the Cou state. Instead, we propose that 
E26 is the major protonation site in the Cou state, because 
E26 is likely to have a higher pK,°“ than D34 according 
to our structural model of the Cow state. Our hypothesis 
is also consistent with the observation that D34 does not 
react with DCCD (a chemical that reacts with proton- 
ated carboxyls) in the Cow state (i.e., in the presence of 
TPP’) [13]. (ii) In the previous model, substrate loading 
occurs after deprotonation in the C,, state. In contrast, 
we propose that substrate loading disturbs the micro-en- 
vironment of D34 and induces its deprotonation. Since 
it is the deprotonation, but not the ligand binding per 
se, that triggers the C,,-to-Co,, conformational change, 
deprotonation would unlikely occur before substrate 
loading. This notion is also supported by the observation 
that D34 reacts with DCCD in the absence of TPP’, thus 
in the Cp state [13] (for more discussion see Supplemen- 
tary information, Data S2). (iii) Most importantly, in our 
model, we explicitly consider the electrostatic interaction 
between the protonation site and the membrane potential 
as the major driving force for the conformational change 
of MdfA (and PMF-driven MFS transporters in general). 

In summary, here we determined the 3D structure of 
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the multidrug-resistance antiporter ecMdfA in complex 
with a known substrate or substrate analogues. The 
structural information combined with functional anal- 
yses improved our understanding of the mechanisms 
of multidrug expulsion by MFS antiporters, including 
substrate-binding modes, competition between substrate 
binding and protonation and functional roles of con- 
served motifs in the transport cycle as well as in ener- 
gy-coupling in general. These new structural information 
may open new avenues to developing novel strategies of 
counteracting drug resistance. 


Materials and Methods 


Expression and purification of recombinant proteins 

The full-length Mdf4 gene was cloned from the genome of £E. 
coli BL21(DE3) and was found to be highly expressed in a GFP- 
based high-throughput screen [41]. An incidental point mutation, 
Q131R, was introduced during the cloning. Since it did not show 
effects in our drug-resistance assay, this mutant was considered as 
WT hereinafter and used in subsequent structural and functional 
studies. MdfA WT and mutant strains were subcloned into a pET- 
28a vector (Novagen) using restriction endonucleases Ncol and 
Xhol. All variant sequences were confirmed by DNA sequencing. 

Recombinant C-His, tag proteins were overexpressed in the 
E. coli C43 (DE3) strain and induced using 0.5 mM isopropyl 
B-d-thiogalactoside (IPTG) at ODsoo nm Of 0.8. After growth at 16 
°C for 18 h, cells were harvested and resuspended in buffer A (20 
mM Tris-HCl, pH 8.0, 300 mM NaCl, 10% (v/v) glycerol, and 5 
mM -mercaptoethanol) and subjected to two runs of homogeniza- 
tion at 10 000-15 000 p.s.i. using a JN-R2C homogenizer (JNBio, 
China). Cell debris was removed by centrifugation at 17 000x g 
for 15 min, and the supernatant was ultra-centrifuged at 100 000x 
g for 1 h. The membrane fraction was solubilized in buffer A sup- 
plemented with 0.5% (w/v) n-decyl-B-D-maltopyranoside (DM; 
Anatrace) for 2 h at 4 °C. After second ultra-centrifugation at 100 
000x g for 30 min, the supernatant was loaded on 2 ml of Ni™”-ni- 
trilotriacetate affinity resin (Ni-NTA; Qiagen) and washed with 
30 ml buffer A containing 20 mM imidazole and 0.2% (w/v) DM. 
The protein sample was eluted with 15 ml buffer A containing 350 
mM imidazole and 0.2% (w/v) DM, and concentrated to about 10- 
15 mg/ml. The concentrated sample was then loaded onto a Su- 
perdex-200 10/30 column (GE Healthcare) pre-equilibrated with 
buffer B (20 mM Tris-HCl, pH 8.0, 100 mM NaCl, 5 mM f-mer- 
captoethanol, 1.2 mM sodium Dxc as well as mixed detergents of 
0.2% n-nonyl-B-D-glucopyranoside (Anatrace) and 0.025% LDAO 
(Anatrace)). Peak fractions were collected, and the pooled protein 
sample was concentrated to 15 mg/ml before crystallization trials. 


Crystallization of MdfA 

Crystals of MdfA were grown at 16 °C using the hanging-drop 
vapour-diffusion method (1 ul + 1 ul over 200 pl). With a reser- 
voir solution (50 M sodium acetate, pH 5.8, 20%-24% (v/v) PEG 
400, 10 mM praseodymium acetate, and 50 mM magnesium ace- 
tate), crystals grew typically to 20 x 20 x 50 um’ in about 3 days. 
Se-Met derivative crystals were grown under the same condition 
as the native ones. Cm-containing crystals were obtained by soak- 
ing the original crystals in the reservoir solution (pH 5.8) supple- 


www.cell-research.com | Cell Research 


Jie Heng et al. hpg 


1071 


mented with 5.0 mM Cm. LDAO-complex crystals were obtained 
by soaking the original Dxc-containing crystal in the reservoir 
solution (substituting sodium acetate with Tris-HCl, pH 8.5). Note 
that LDAO in the protein sample was at 0.025% (~1 mM) concen- 
tration. Crystals were flash-cooled in liquid nitrogen for storage 
and data collection. 


Data collection and structure determination 

X-ray diffraction data sets were collected at KEK or SSRF and 
processed with the HKL2000 package [42]. WT MdfA contains 
18 Met residues. Positions of sixteen of the seleno-substituted 
Met (SeMet) residues were identified in the initial anomalous dif- 
ference Fourier map of the MdfA-Dxc complex crystal and were 
used in phase determination by using the single-wavelength anom- 
alous dispersion (SAD) methods implemented in the program 
HKL2MAP [43]. The experimental SAD phases were of excellent 
quality so that ~80% of the structural model was built automatical- 
ly by the computer program Phenix.autobuild. The space group of 
the crystal form is C2, and there is one MdfA molecule per crys- 
tallographic asymmetric unit. The model was further refined using 
the program Coot [44]. Model validation was carried out using the 
web-based program Molprobity [45]. 


Drug-resistance assays 

E. coli C43 (DE3) [46] strain was transformed with the MdfA 
gene-containing pET28a plasmids, individually. Single clone was 
picked from LB plates to 48-wells plate containing 1 ml Overnight 
Express Instant TB Medium (Novagen). After 5-h propagating at 
37 °C, cultures were transferred to 16 °C and grown further for 10 
h. Cells were harvested, resuspended and adjusted to OD, of 2.0, 
sequentially diluted as specified, and spotted over LB agar plates. 
The growth medium was supplemented with 30 ug/ml kanamycin 
(selection marker), 0.2 mM IPTG and 2.5 ug/ml Cm (resistance 
marker). The ability of the E. coli cells to form colonies was re- 
corded by imaging of the plates after 37 °C incubation for 36 h. 
The remaining cells were sonicated for 3 min and centrifuged at 
17 000x g for 10 min. The supernatants (5 u1), which contained 
the membrane fraction, were subjected to (10%) SDS-PAGE. An- 
ti-glyceraldehyde 3-phosphate dehydrogenase mouse monoclonal 
antibody and anti-His mouse monoclonal antibody (Transgen Bio- 
tech) were used to detect the expression levels of GAPDH (internal 
control) and MdfA variants. 


ITC assays 

The binding affinities between ligands and MdfA variants were 
measured with an ITC200 micro-calorimeter (MicroCal). WT and 
mutatant variants of MdfA were prepared in a buffer containing 
40 mM HEPES (pH 7.2), 100 mM choline chloride, 0.01% (w/ 
v) lauryl maltose neopentyl glycol and 0.002% (w/v) cholesteryl 
hemisuccinate Tris salt. Cm was dissolved in the same buffer. Pro- 
tein concentration in the measurement cell was 50 uM, while the 
ligand concentration in the titration procedure were adjusted to 2-5 
mM. Titrations were performed at 25 °C, and data were fitted to 
the one-binding-site model using the software Origin 7.0 (Micro- 
Cal). 


Protein Data Bank accession number 

The diffraction data and refined coordinates of the crystal struc- 
tures of MdfA are deposited into PDB. The accession codes are 
4ZOW for the MdfA-Cm complex, 4ZP0 for the MdfA-Dxc com- 
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plex and 4ZP2 for the MdfA-LDAO complex. 
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